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Abstract. We propose that a weakly-coupled nonlinear dielectric waveguide - 
surface-plasmon system can be formulated as a new type of Josephson junction. Such a 
system can be realized along a metal - dielectric interface where the dielectric medium 
hosts a nonlinear waveguide (e.g. fiber) for soliton propagation. We demonstrate 
that the system is in close analogy to the bosonic Josephson- Junction (BJJ) of atomic 
condensates at very low temperatures, yet exhibits different dynamical features. In 
particular, the inherently dynamic coupling parameter between soliton and surface- 
plasmon generates self-trapped oscillatory states at nonzero fractional populations with 
zero and tt time averaged phase difference. The salient features of the dynamics are 
presented in the phase space. 

PACS numbers: 05.45.-a, 42.65.-k, 03.75.Lm 
1. Introduction 

Plasmonics encompasses the science and technology of plasmons which are the collective 
oscillations of electrons mainly in metals. In particular, the coupling of plasmons with 
photons give rise to a hybrid quasiparticle known as the surface plasmon-polariton (also 
called surface-plasmon) pj, which can propagate along metal surfaces. As a result, light 
can be coupled to, and propagate through structures that are much smaller than it's 
wavelength. Thus, plasmonics enables the subwavelength photonics [21 13] and bridges 
it to electronics at nanoscale and offers promising applications in nano-optics and elec- 
tronics such as lasing, sensing [H [5l [6l [7] . 



In this context, the coupling between the surface plasmon polaritons and the 
dielectric waveguide modes is of interest for integrated optoelectronic applications. 
Much effort has been devoted to understand and enhance this coupling in guided mode 
geometries |8l HI [10] and in layered systems [11] . In particular, a recent work describes 
the resonant interaction model between the surface-plasmons on a metal surface and a 
soliton in a nonlinear dielectric medium. [12] In that model, the interaction depends 
on, and hence can be controlled by the soliton amphtude which leads to coupled 



The nonlinear dielectric waveguide - surface-plasmon Josephson Junction 



2 



surface-plasmon-soliton modes, and provides a novel way of manipulating the surface 
plasmon propagation. Motivated by this study, we propose in this paper that a weakly 
coupled nonlinear dielectric waveguide surface-plasmon system can be formulated as 
a new type of 'dielectric waveguide - surface-plasmon Josephson Junction (DWSP- 
JJ)' and investigate its dynamics within the model introduced in Ref. [12] • This 
formulation enables a remarkable connection between this optical system and the well- 
known Josephson- Junction dynamics of Bose-Einstein atomic condensates [13] and opens 
a new perspective to quantum plasmonics, where novel quantum optical phenomena can 
be reahzed and utilized via coupling to surface plasmons. [H] 

The paper is organized as follows. In Section [21 the theoretical model is introduced. 
Section [3l investigates the dynamical features of the DWSP-JJ system in comparison to 
SJJ and BJJ systems. Section H] concludes the paper. 



2. Model and Theoretical Formulation 



Within the scope of this paper, we adopt the model system given in Ref. [12]. This 
intuitive model consists of the co-propagating optical soliton and surface-plasmon (SP) 
electric fields along a metal-dielectric interface [201 |2T]. The nonlinearity in the dielectric 
is assumed to be confined at a distance, d, from the metal interface such that the surface- 
plasmon propagation retains linearity, and the weak coupling between the soliton and 
SP-fields can be treated perturbatively. Assuming translational invariance in the z 
direction, and choosing y as the propagation direction, the total electric field of the 
coupled system is written in the form, 

E(..,) = c,fe)e-- + ^^j;^^^. (1) 

where Cp^siy) are the surface-plasmon and soliton amplitudes, multiplied by the 
respective transverse profiles with = ^k^ — k"^ and = fc-y/7/2|cs|. The propagation 
wavevector is k whereas kp is the surface-plasmon wavevector. The nonlinearity 
parameter of the dielectric medium is 7. This product ansatz for the lateral and 
longitudinal field profiles reduces the system to essentially a one-dimensional problem, 
where the amplitudes Cp^s{y) obey the following coupled oscillator equations, that are 
written in the dimensionless coordinate ( = ky 

Cp + f3^Cp = q{\cs\)cs, 

Cs + (3!cs = q{\cs\)cp. ^ ^ 

Here, f3p = kp/k and f3s = 1 + 7|csp/4 are the propagation constants. The coupling 



parameter is given by 5'(|cs|) ~ exp{—ky 'y /2\cs\d) [12j. The rationale of this functional 
form is that the lateral tail of the soliton field, which depends on the soliton amplitude 
\cs\, acts as a source to excite the surface plasmons. Substituting Cp^s = Cp,sG^'' and 
employing the slowly varying amplitude approximation, the following set of equations 
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are obtained for the amplitudes Cp^s 

_jn — n _ <i{\CB\) fi 

= -^Cp + z/,(|C,|)C,. ^ ' 

In Eq. [Sli/p^/Sp — 1<^1, i/s = /3s — 1^1 are the small deviations of the 
dimensionless propagation constants of surface-plasmon and soliton, respectively. The 
eigenmode analysis indicates that the resonant coupling g^es = exp[— /c(i^2i/pj occurs 
around z/p = 0.2, 7] = 0.2 [12]. 



2.1. Josephson Junction Formulation 

Remarkably, the Eq. [3] can be cast in analogy to the Josephson junction dynamics 
by writing Cg^p = C^^pC^'^^-'', and introducing the fractional population imbalance 
Z = (|Csp — \Cp\'^)/N and the relative phase difference between the soliton and the 
surface plasmon = 0^ — 0p as follows: 



Z = -q{Z) Vl - Z^ sin 0, (4) 

= + + 4^£1^ cos 0. (5) 

We set = (|Csp + |Cpp) = 1, a normalized constant for the isolated system with 
no population dissipation, and define the parameters A = ri/2., AE = ri/2 — Up, with 
rj = A is the nonlinearity (i.e. the soliton strength) whereas AE parametrizes the 
asymmetry between the soliton and surface plasmon states occupied by the photons. The 
coupling parameter takes the Z-dependent form 



q{Z) = 6-^=^2^(1+^). (6) 

An immediate comparison reveals the similarities and differences to the dynamical 
model of the bosonic Josephson junction (BJJ) in a double well trap model of two 
Bose-Einstein condensates 



Z = ~Vl - Z2sin0, (7) 

= AZ + AE + -^=£=cos0 (8) 

In the BJJ model, A describes the interatomic interactions and AE is the difference 
between the zero-point energies of the trapping wells. Both A and AE are dimensionless 
parameters scaled by the coupling matrix element. [22] We observe that the DWSP-JJ 
equations would reduce to that of the BJJ model for a constant coupling parameter 
q{Z) = 1 provided that the other parameters are scaled accordingly. 

It is instructive to present also the case of superconducting Josephson junctions 
(SJJ) where the dynamical variable is the voltage across the resistive shunt between the 
two superconductors [T8] : 

(9) 
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The SJJ tunneling AC current is given by, 
/ = Jo (cos 0)0. 



(10) 



The Cooper-pair population imbalance across SJJ is zero when the SJJ is closed over 
an external circuit [17]. In the absence of external circuit, an isolated SJJ can exhibit 
coherent Cooper-pair oscillations only at very small amplitudes [171123]. By virtue of the 
oscillating current, the SJJ is generally discussed in terms of a rigid pendulum analogy. 
On the other hand, the BJJ dynamics resembles to that of a nonrigid pendulum with 
length-dependent angular momentum (Z) [22]. 

In the BJJ, the double-well trap is created by a laser field barrier that divides a 
single trapped condensate in two parts. Hence, the asymmetry of the wells {AE), as 
well as the barrier height (i.e. coupling) is controlled externally by the laser field. The 
coupling can further be modulated temporally by the laser beam. 

To this end, we stress that the presence of Z-dependent nonlinear coupling in 
DWSP-JJ implies fundamentally different dynamical features from the BJJ and SJJ 
models, particularly in certain parameter regimes, as discussed in the next section. 

3. Results 

We start exploring the dynamical landscape of the DWSP-JJ model first by determining 
the stationary solutions at fixed points. Equation (jl]) shows that the fixed points occur at 
(0*=27™) or (0*=(2r;,+l)7r) with the respective population imbalance values determined 
from 



We investigate the range 0.01 < Up < 0.35, and 0.01 < rj < 0.35 including the 
resonant coupling [vp = 0.2, rj = 0.2) [12]. For the scaled distance parameter kd we use 
the range 3-12 and discuss the uncoupled limit {kd » 1 qualitatively. 

3.1. Fixed Points 

The zero-phase modes describe the transfer of energy between the soliton and surface- 
plasmon states with zero time-average value of the phase. In figures [Tj^a-c) the f(Z) is 
plotted for different values of AE for Vp = 0.15 and kd = 3, 6, 12. For AE = the fixed 
point is at Z = 0, which can also be deduced by inspecting Eq. ( ITT]) . This is analogous 
to the "symmetric double well" of the BJJ model. For AE < (>)0, the fixed point 
occurs at nonzero fractional population imbalance Z > (<)0. 

The other mode is the vr-phase mode, where the time-average value of phase is 
< >= vr. Up to three fixed points emerge depending on the model parameters as 
shown in Fig. The location of the fixed points depend strongly on the scaled dis- 
tance kd. 




q{Z*)Z* 




(11) 
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Figure 1: f{Z) (Eq. [TTJ for zero-phase modes plotted at Up = 0.15 with various values of /\E and (a) kd = 3 (b) 
kd = 6 and (c) kd = 12 The critical points are the roots of f{Z) = 0. 




Figure 2: f(Z) (Eq. [TTJ for 7r-phasc modes plotted at i/p = 0.15 with various values of Ai? and (a) kd = 3 (b) kd = 6 
and (c) kd = 9 The critical points arc the root(s) of f{Z) = 0. 

3.2. Phase-space analysis 

A concise description of the dynamical behavior can be given by the phase-space rep- 
resentation of Eq. lU In Figure [31^ a-c), several phase space trajectories are given for the 
symmetric case with z/p = 0.15 and kd = 3, 6, 12. For kd = 3 (Fig. [3]^a)), the zero-phase 
and TT-phase modes exhibit similar behavior, with rigid-pendulum type closed-orbit oscil- 
lations at small amplitudes and anharmonic closed-orbit oscillations at large amplitudes. 
The running phase trajectories consists of 1^1 ~ 1 plateau connected by a population 
inversion occuring close to ~ {2n + 1)71/2 points. All these features are similar to the 
dynamical modes of the BJJ model discussed elsewhere [T5] . 

When the spacing parameter is increased to kd = 6, we observe a drastic change 
in the phase-space portrait (Fig|3t^b)). The anharmonicity of the large-amplitude zero- 
phase modes become prominent. For the vr-phase, three fixed points appear, of which 
the Z* ~ 0.98 and Z* ~ —0.25 are enclosed by bounded trajectories, whereas the 
Z* = is an unstable fixed point. In Figure [3]^c), we observe that increasing the kd 
further {kd =12) results in the decoupling of the system for soliton-dominant initial 
populations {Z{0) > 0.5). This is expected since the coupling parameter q{Z) remains 
small. On the other hand, a surface-plasmon dominant initial population {Z{0) < 0) 
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is still effectively "strong coupled" since the coupling parameter can be close to unity 
even though kd is large (see Eq. (|6])). The closed-orbits of the zero phase mode widens 
in phase towards the two unstable fixed points as opposed to being confined around 
1^1 < 7r/2. The closed orbit vr-phase modes occurs around Z* ~ —0.8. In accordance 
with Fig. W^c) another fixed point exists that occurs almost at Z* = 1 and invisible in 
this scale. 




Figure 3: Several phase-space trajectories of DWSP-JJ system for AE = 0,Up = 0.15, (a)fcd = 3, (b) kd = 6, (c) 
kd = 12. 




The phase trajectories for {AE ^ 0) and fixed Up = 0.15 are given in Fig. HJ 
Figure H] (a) is plotted for AE = —0.025, kd = 6. The zero-phase fixed point is 
located at Z* = 0.19. The vr-phase closed modes occur around Z* = —0.6 and in a 
narrow region around Z* ~ 1. For 0(0) = running phase orbits have large amplitude 
oscillations for 0.81 < Z < 0.85 and small oscillations for 0.85 < Z. The phase diagram 
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for AE = 0.025, kd = 6 is shown in FigH] (b). The zero phase mode is at {Z* = —0.2). 
The TT-phase mode has closed orbits only around {Z* = 0.95) . When kd is decreased 
{kd = 3) (FigJl](c)) = vr mode forms closed orbits around Z* = 0.9. 

In the BJJ model, the parameter AE is a measure of the asymmetry between the 
two trapping states of atoms. Evidently, a nonzero AE induces a population imbalance 
between the two states. In the DWSP-JJ model, this asymmetry is indicated by the 
location of the (p* = 0, Z* fixed points. 

We next investigate the variation of Z{Q, q{Z), and 0(C) for various trajectories 
from Fig ll](a). Figure [5] show the propagation with initial relative phase 0(0) = 
for different initial values of Z{0). As expected from the functional form in Eq. ([6]), 
the coupling parameter shows variations commensurate with that of the population 
imbalance. Small amplitude oscillations (Fig. [5]^a), circles and dotted line) have small, 
almost constant q{Z) whereas large amplitude oscillations (triangle and square) exhibit 
impulsive behavior of the coupling. Similar observations are present for 0(0) = vr in[6l 
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We conclude our discussion by a qualitative comparison of the BJJ, DWSP-JJ and 
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SJJ phase portraits in Fig. [7] (a-c), respectively. The normahzed SJJ AC current is 
plotted for a constant junction potential in Eq. |9l A detailed comparison between SJJ 
and BJJ was reported in the literature [T5]. Here we rather focus on the BJJ and 
DWSP-JJ phase space features. 

The BJJ model is depicted with an effective coupling constant of g ~ 0.03 which is 
within the dynamic range of the q{Z) of DWSP-JJ. The phase space of BJJ is decorated 
with open- and bounded-phase trajectories that are symmetric with respect to the Z 
axis. The 0(0) = trajectories are bounded for |2'(0)| < 0.8. The bounded orbits 
are harmonic for |2'(0)| < 0.5 and become anharmonic for 0.5 < Z{0) < 0.8. The 
0(0) = vr trajectories are open for |2'(0)| < 0.9 and become bounded anharmonic orbits 
for 0.9 < \Z{0)\ < 1. 

In the DWSP-JJ phase portrait plotted with a comparable parameter set 
(Fig. [71(b)), the asymmetry in Z induced by the dynamical coupling parameter is 
prominent. Qualitatively speaking, the DWSP-JJ acts like a double-well trap BJJ 
system, where the barrier height (i.e. the coupling) between the two wells depends 
on the population of one well (i.e. the soliton amplitude). The bound trajectories are 
strongly anharmonic except for \Z{0) < 0.2 1. For 0(0) = 0, the orbits are closed for all 
negative Z{0) values and for positive Z{0) < 0.55. The 0(0) = vr states are bounded for 
-0.9 < Z{0) < 0. 




Figure 7: Comparison of the phase-space trajectories of (a) BJJ (q = 0.03, AB = 0, A = 0.15) and (b) DWSP-JJ 
{kd = 6, AE = 0, A = 0.15) and (c) SJJ models 

Finally, we would like to point out the dissipative effects that are excluded from the 
present dynamical model. Naturally, both population- and phase-dependent dissipation 
mechanisms are applicable to this model. The effects of phase dissipation are discussed 
in part in [12]. We plan to discuss a dissipative- DWSP-JJ model in a subsequent work. 
In particular, the absorption loss of the surface-plasmon will be relevant and may alter 
the phase portrait substantially for strong loss. Nevertheless, the qualitative features of 
the present analysis are still providing insight and reveal intriguing properties that can 
be observable under weak dissipative conditions. 
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4. Conclusion 

In this work, we proposed a weakly coupled optical-soliton and metal surface-plasmon 
system as a novel type of Josephson Junction and investigated its dynamical properties. 
We have found that the coupling parameter that depends on the soliton amplitude allows 
rich dynamical features, different than that observed in bosonic Josephson Junctions. 
The DWSP-JJ is potentially a convenient system to investigate collective dynamics of 
photons and surface plasmons, that may reveal novel classical and quantum plasmonic 
phenomena. 
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